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Functional relations for the density functional exchange and correlation 
functionals connecting functionals at three densities. 

Daniel P. Joubert 1 '^ 

Centre for Theoretical Physics, University of the Witwatersrand, PO Wits 2050, Johannesburg, 
South Africa 

(Dated: 21 September 2011) 

It is shown that the DFT exchange and correlation functionals satisfy 

= jE hx [ PN ] + 2E2 [ P n] 

-lE hx [pl_ 1 ]-2E2[pl_ 1 ] 

+2 / dV (p^i (r) - Pl- X (r)) v° ([p N ] ; r) 



+ / d 3 r' (p ^ (r) - pl_, (r)) r.V^ (\p N ] ; r) 
ofr'pN (r)r.V^([ PJ v];r) 

d 3 r'P (r)r.V^ c (M;r) 



rf 3 ^/ 7 (rK xc (M;r) 

In the derivation of this equation the adiabatic connection formulation is used where the ground state density 
of an iV-electron system, pn, is kept constant independent of the electron-electron coupling strength 7. Here 
Ehx [p] is the Hartree plus exchange energy, EJ [p] is the correlation energy, v^ xc [p] is the Hartree plus 
exchange-correlation potential, v c [p] is the correlation potential and v [p] is the Kohn-Sham potential. The 
charge densities pm and pjj_ 1 are the N- and (N — l)-electron ground state densities of the same Hamiltonian 
at electron-electron coupling strength 7. p (r) = px (r) — P^-x ( r ) is the Fukui function. This equation can 
be useful in testing the internal self-consistency of approximations to the exchange and correlation functionals. 
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Density Functional Theory (DFT)i is one of the most 
important tools for the calculation of electronic struc- 
ture and structural properties of molecules and solids. 
In all practical applications of DFT, however, approxima- 
tions to the exact functionals have to be used 2-6 . Exact 
relations for density functionals and density functional 
derivatives can play an important role in the development 
of accurate approximations to the exact functionals. A 
successful approach to the design of improved approxi- 
mate density functionals is by 'constraint satisfaction' 7 , 
where the approximate functionals are required to sat- 
isfy properties of the exact functionals. The accuracy 
of approximate functionals can be tested by comparing 
to accurate calculations or to experimental data. A use- 
ful additional test will be to determine the internal self- 
consistency of approximations. With this in mind, the 



following expression is derived: 

= jE hx [p N ] + 2E2 [p N ] 

- lEhx [pl_ 1 ]-2E2[pl_ 1 ] 

+2 I d 3 r> (p%_, (r) - pl_ x (r)) v° ([p N ] ; r) 



+ / d 3 r' (pk-i (r) - P' N ^ (r)) r. W Qp N ] ; r) 
ofr'p N (r)r.V^([ PJV ];r) 
dh'pl^^r.VvJdpl^ir) 
d 3 r'P (r)r.Vvl xc ([p N ];r) 
2 I dV/ 7 Wi c (M;r) 
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In the derivation of this expression the adiabatic con- 
nection formulation 8-11 is used where the ground state 
density of an ./V-electron system, p^, is kept constant 



2 



independent of the electron-electron coupling strength 
7. In Eq.fll]) v 1 hxc [p] is the Hartree plus exchange- 
correlation potential, v c [p] is the correlation potential 
and v° [p] is the Kohn-Sham potential. EJ [p] is the cor- 
relation energy while Eh x [p] is the Hartree plus exchange 
energy. The charge densities pn and p 1 N _ 1 are the N- 
and (N — l)-electron ground state densities of the same 
Hamiltonian at electron-electron coupling strength 7 as 
discussed below, p (r) = p^ (r) — p~' N _ 1 (r) is the Fukui 
function. 

Equation (TTJ) couples functionals and functional deriva- 
tives evaluated at three different densities and two par- 
ticle numbers. This equation can be used as a stringent 
test to check the internal self-consistency of approxima- 
tions to the exchange and correlation functionals. If exact 
exchange is used it can serve as a check on approxima- 
tions to the correlation energy functional. 

Evaluation of ([TJ) requires two independent self- 
consistent Kohn-Sham calculations, one for the N- 
electron system to determine pn and p%_i, and one for 
the (N — l)-electron system to determine pjf_ v Once 
the densities are found, the functionals and functional 
derivatives can be used in Eq.([T]) as a test of how well 
the approximations work. 



I. PROOF OF EQUATION 

According to the Hohenberg-Kohn theorem 1 - and its 
generalization to degenerate states 12 the ground state 
energy E 1 [p] of a system of interacting electrons is a 
functional of p, the ground state densit y 13 ' 1 . 

BP [p] = T 1 [p] + 7KI [P\ + fd 3 rp (r) ( [p] ; r) (2) 



where T" 1 [p] is the kinetic energy and jV^ e [p] the mu- 
tual Coulomb interaction energy at density p. In order 
to derive Eq.([T]) the adiabatic connection approach will 
be used 8 - - — in which the external potential uj xt ([p] ; r) 
is constructed to keep the ground state density indepen- 
dent of the electron-electron interaction strength, scaled 
by 7, and has the for m 15 ' 16 



= (l-j)v hx {[p];r) 

+ «c(M;r) -v<{[p];r)+v cxt {r). 



(3) 



u 1 ([p] ; r) = v cxt (r) is the external potential at full cou- 
pling strength, 7 = 1, and u ([p];r) is non- interacting 
Kohn-Sham potential. The exchange plus Hartree 
potentia l 13 ' 14 Vh x ([p];r), is independent of 7, while the 
correlation potential ([/"]; r ) depends in the scaling pa- 
rameter 7. The adiabatically scaled TV-electron Hamilto- 
nian H 1 has the form 8-11 



Atomic units, % = e = m = 1 are used throughout. T is 
the kinetic energy operator, 



1 N 



and "/V oc is a scaled electron-electron interaction, 

N 

iVee = 7 ' 



1 



and 



i<3 



N 



v 1 [p] = [p];ri. 

1=1 

Define the energy functiona l 13 ' 14 

F 1 \p\=Ti [P]+7K1[P] 

= T° [p] +1 E hx [ P ]+E2[p] 



(5) 



(6) 



(7) 



(8) 



where Eh x [p] is the Hartree plus exchange energy, and 
El [p] is the correlation energy. Note that Eh x [p] is in- 
dependent of 7. As shown in appendix |A"1 



= »•» (p) + -r.V.l (r) + - («J, t (r) - ij (p)) 



2 J PV Sp(r')Sp(r) 



(9) 



where tj ([p] ; r) = & Jp^ with Xj? [p] the correlation part 

of the kinetic energy and v7 (r) = -fVh x (r) + (r) . 
For notational convenience the functional dependence on 
p has been suppressed. The chemical potential p depends 
on the asymptotic decay of the charge densit y 17 ' 18 , and 
hence, since the charge density p is independent of 7 by 
construction, p is independent of 7. 
The Fukui function^^ 



r M- 



5p 



5v~< (r) 



Sp(r) 



N 



SN 



satisfies 



dp 
dN 



6p(r)6p(r>y (r } 
J d 3 r'r (r') = 1. 

J d\' (3p(r')+r".Vp(r')) = 0, 
it follows from ©, (QT]) and (pj} that 



and 



Since 



(10) 



(11) 



(12) 



(13) 



p = / d 3 r'P (r) (r) + - / dV P (r) r.Vu 7 (r) + 



[p\- 



(4) 



- / A7 7 (r)(t(r)-?(r)) 



(14) 
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Using ([3]), this equation can be reduced to 

d 3 r'P (r) [2v° (r) + r. W (r)] - 2p 

d\'P (r)r.Vt^ c (r) 

f / d 3 rT (r) K K W+t?(r)]. 



(15) 



From the virial theorem for the Kohn-Sham independent 
particle wavefunctions^ 1 - 



dV [ PN (r) - p^-i (r)] [2u° (r) + r.Vv (r)] 
= 2(£°M- J E [p^_ 1 ]), (16) 

where £° [p] and £° [p^_ x ] are the Kohn-Sham ground 
state energies of the independent N- and (N — l)-particle 
systems of the same Kohn-Sham Hamiltonian with po- 
tential v° . By definition, since p is independent of 7 17 ' 18 , 



M = EP [p] EP K_J 



(17) 



where E 1 [p N ] and E 1 [p]^_ 1 ] are the ground state en- 
ergies of the N- and (N — l)-particle systems of the 
same Hamiltonian i? 7 , Eq. Q, with p^ and the 
corresponding ground state densities. By construction 
P~n = P N i s independent of 7 but pjj_ 1 is expected to be 
a function of 7. 

Since the Fukui function 22,23 



P (r) 



Mr) 



<S7V 



Pn 0) - Pat-i (r) 



(18) 



is the difference between the ground state densities of the 
N- particle and (N — l)-particle systems, it follows from 
(fl5l). (fl"6]) and El, that 



2 yd 3 / (p^.iW-p^W)* (r) + 
+ I dV ^^(rj-^fr)) r.W (r) 

dV/ 7 (r)r.V^ c (r) 
- / rf 3 r'f(r)K IC (r)+?W] (19) 



The reference to the correlation part of the kinetic energy 
can be eliminated as follows. Sincei 5 - 



7K 7 e [P] = lE hx [p] + EJ [p] - T-' [p] 
the last line in (IT9"]) can be written as 

d 3 r'P W K xc (r) + t2 (r)] 
+2 fd s rT(v)vl xc (v) 



(20) 



As shown in appendix [B? 4 . 



dV/7 (r) ijrpr = vg [pn] ~ vg [pl - i] ■ (22) 



,15 



Sinci 



£ c 7 [P\ + T2 [ P } = - d'r'p (r) r.V [vj ([p] ; r)] , (23) 



it follows from j2Q]), (j2Tj), d22J) and ([23]) that 

d 3 r'/ 7 (r) bL c W + C (r)] 
-! / d 3 r'p N (r)r.Vv-<([pN};r) 



7^ K-i] + 2i? c 7 K_i] 
-lE hx [p N ] - 2£ 7 [Pw] 

+ 2 / ATWieW- 



(24) 



Combining (19]) and ([24]) leads to Equation (T]). 

This equation is now entirely in terms of the exchange 
and correlation energy and potentials and other known 
density dependent quantities. 



II. DISCUSSION AND SUMMARY 

Equation (1]) is valid for pure state DF T 13 i 14 since equa- 
tions (Til?]) and ([2"2")) have been shown to be correct for pure 
states. Throughout it was assumed that the functionals 
derivatives are well define d 13 ' 14 ' 25 . 

Note that p]q_i ^ P%-i- In order to determine pTv-i 
a Kohn-Sham calculation with 26 

V ° [PN-l] = V ~' [PN] + IVhx Kr_i] + v~c [Pn-i] 

has to be performed. Note that v 1 [pj 4 _ 1 ] = v 1 [pn] ■ 

In summary, an equation that couples exchange and 
correlation functionals and functional derivatives evalu- 
ated at three different densities and for two particle num- 
bers has been derived. This equation can be used as a 
stringent test to check the internal self-consistency of ap- 
proximations to the exchange and correlation functionals. 
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Appendix A: Derivation of Eq. J9) 

Let p\ (r) = X 3 p (Ar) , the uniformly scaled density. 
Then (see for example Eq.(A.33) in Appendix A of 
reference^ 3 .) , 



d SF^[ Px ] 



dX Sp x {r) 



Pa,A=1 



dV (3p (r')+r'.V'p (r')) 



5 2 F'[p] 
8p{v')5p{r) 



= -jdh'p{v')v'.V * 2F ' [,l] 



Sp (r') Sp (r) 



(Al) 



The last line is valid if p (r') vanishes when r ~ > oo as 
would be the case for a finite system. Now consider the 
Schrodinger equation 



JV JV 



N 



J2 yx (W; r *) 



A f—f |r, - r., 



({*}) 



from which it follows that 

JV 



(«* [p]) *x ({ r .}) 
ows that 

JV JV 



(A2) 



(=i 



JV 



i=l 

A 2 £* 



*x ({ArJ) 



(wi *t ({Ar,}) 



(A3) 



If normalised ^ * ({r^}) yields p (r) then nor- 
malised ^ a ({Ari}) yields A 3 p(Ar). The Levy 
constrained minimization approach 12 implies 
that normalised ^~ ({Ar-j}) yields A 3 p (Ar) 



T + -iVe, 



and minimizes ( ^ 

' *^A 3 p(Ar) 

Thereforei^^T {[p x ] ; r) = X 2 vi ([p] ; Ar) and 



S P x (r) 



+ X 2 v~ {[p];Xr)=fil, 



(A4) 



and 



dX 



(IPA];r) 



A=l 



= 2v-> ([p] ; r) + r.W ([p] ; r) - ^ Qp] ; r) (A5) 
For an iV-electron system, 

^l,x = El (X 2 vi [p]) - El_, [X 2 vi [p]) (A6) 



with Ejf 



^X 2 vl [p N fj and El_ x \X 2 vi [p N ]j the 
groundstate energies of the N- and (N — 1)- electron 
systems with the same external potential A 2 t>x [p] . From 
Eq.(lA"3l it follows that 



Mjv.a = El (X 2 vi - El_, (X 2 vi [p]) 

= X 2 (El(vi [p])-I<$-i (yi [p])) 

= X 2 p (A7) 

where the last step follows from the definition of the 
chemical potential of the jV-electron system. Note that 
the chemical potential p is independent of 7 17 ' 18 . 
The correlation energy E% [p] is defined asi 3 - - — 



El [p] 



p 



T + >yV e , 



T + 7 ^e< 



(A8) 



where is the ground state wavefunction of H 1 

that yields p and \^° p ) the Kohn-Sham independent JV- 
electron groundstate wavefunction that yields the same 
density. The correlation part of the kinetic energy is given 
by 



T 



T 



q, 



and "fV^, [p}^ — can be written as 

jV e l [p] = 1 E hx [p] + E2 [p] - T2 [p] 



(A9) 



(A10) 



where E/- IX [p] is the sum of the Hartree and exchange 
energies. Note that Eh x [p] is independent of 7. It follows 
from (EH1) and JM) that 



(All) 



tl (r)+l- r v2(r)-v2(r)=0 } 
a 7 



where tj ([p] ; r) = and v< ([p] ; r) = C 



J5p(r 



bining ©, ([AT]). ([A~4) . (1A31) . (IA71) and (lATTj) yields Eq. 



Appendix B: Derivation of Eq. {[22) 

Multiply Eq.© by the Fukui function and integrate 
over r. Use Eq. fTTj) , the fact that 



dV(3p(r') + r.Vp(r')) = (Bl) 

and Eq. (|All[ ) to arrive at 

dV [2v^ ({p} ; r') + r.W ([p] ; r)] P (r) + 

d 3 r' [ 1Vhx ([p] ■ r') + v2 ([p] ; r') - % ([p] ; r')] f (r') 
2p. (B2) 
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From the definition of EJ [p] and EJ [p] , Eqs. (JMJ and 
191) and the fact that 



Eh* [ P ] = h° p 



(B3) 



it follows that 

,«7e [P] 



7- 



(5p (r) 



7 ^([p];r)+^([p];r)-^([p];r). (B4) 



Using the virial theorem^ and p = i? 7 [pjy] —E 1 \(P N _^\ , 

dV [2^([p];r')+r.V^([p];r)]r (r) 
= 2 M - 7KI [PJV] + jV e l [p^_J . (B5) 

Combining Eq. ([Bljl. (lB4f and (JBSJ) yields Eq. ([22]). For 

a detailed discussion see reference^. 
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